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A model is studied which describes a propagation of a subpicosecond optical pulse in dispersion- 
managed fibre links. In the limit of weak chromatic dispersion management, the model equation is 
reduced to a perturbed modified NLS equation having a nonlinearity dispersion term. By means of 
the Riemann-Hilbert problem, a perturbation theory for the soliton of the modified NLS equation 
is developed. It is shown in the adiabatic approximation that there exists a unique possibility to 
suppress the perturbation-induced shift of the soliton centre at the cost of proper matching of the 

\Q • soliton width and nonlinearity dispersion parameter. In the next-order approximation, the spectral 

f^*) ' density of the radiation power emitted by a soliton is calculated. 
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I. INTRODUCTION 

\Q . The dispersion management for shortpulse propagation in optical fibres is the key current technology for ultrafast 
' high-bit-rate communication lines 0, 0> 0, 0> S B • The basic idea behind the dispersion management borrowed 
from linear systems |10| consists in a compensation of the chromatic dispersion by means of periodically incorporating 

HH additional fibre sections with opposite sign of chromatic dispersion. Optical pulses in a fibre link with incorporated 
sections, the so-called dispersion- managed solitons (DM-solitons), exhibit a number of true-soliton-like properties, 
including elastic scattering and high stability. Such a behaviour is unexpected because the underlying equation, the 
nonlinear Schrodinger (NLS) equation with variable chromatic dispersion coefficient, is by no means integrable. A 
recent example demonstrating a great success of the DM technology is a development of the commercial DWDM 

l— ~~ communication line in Australia with 160 frequency channels and 1.6 Tbit/s total data transmission rate over 3875 
I , km. 

In treating analytically the DM NLS solitons, asymptotic methods have been proved to be especially effective 
0^ ! H, IE Hj|. Using the multiscale asymptotic expansion, Gabitov and Turitsyn decomposed the pulse dynamics 
■ into rapidly varying phase and slowly evolving amplitude and derived an integro-differential equation (the Gabitov- 
Turitsyn equation) for the pulse amplitude. It was shown p^ . ll3llT3| that in the limit of weak dispersion management 
this equation is reduced to the perturbed NLS equation. 

Up to now, DM solitons in fibres are considered within the NLS-based systems, i.e., for picosecond optical pulses. 
The demand for an increased bit rate of a communication line can be implemented by using pulses of subpicosecond 
duration (^100 fs). For ultrashort pulses, the NLS equation becomes inadequate because more subtle effects, the main 
ones being the nonlinearity dispersion, the Raman self-frequency shift and linear third-order chromatic dispersion, 
come into play. It is important that, even for small strengths of these new effects, some of them cannot in general 
be considered as perturbations of the NLS equation. It is precisely such a situation that arises when attempting to 
^ ■ account for the nonlinearity dispersion. On the one hand, it is well known jlaj that the NLS soliton being subjected to 
the action of the nonlinearity dispersion, exhibits the so-called self-steepening. On the other hand, the NLS equation 
with the additional nonlinearity dispersion term, the modified NLS (MNLS) equation, is still integrable and possesses 
its own solitons that propagate without distortion. The reason for such a difference lies in the fact that the true 
MNLS soliton is non-perturbative with respect to the nonlinearity dispersion parameter and cannot be obtained from 
the NLS equation with the nonlinearity dispersion term as a perturbation. 

In the present paper we give a self-contained exposition of the MNLS soliton dynamics in fibre links with variable 
chromatic dispersion. In section 2 we derive the Gabitov- Turitsyn-like integro-differential equation for the soliton 
amplitude and reduce it to a perturbed MNLS equation in the limit of weak chromatic dispersion management. In 
section 3 we construct the soliton solution of the MNLS equation using the Riemann-Hilbert (RH) problem. The 
MNLS soliton solution has been previously derived by different methods [l6L [l7l Il8j . The representation for the 
MNLS soliton given here is the most simple and transparent. Section 4 is devoted to a development, within the RH 
problem approach, of a perturbation theory for the MNLS soliton. Note that the applicability of the RH problem for 
treating soliton perturbations has been initiated in papers |l9l l2('| . We considerably simplify the formalism and final 
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expressions, as compared with the previous results |2lLl22ll23l | on the perturbed MNLS soliton. In section 5 we work 
out in detail the general results of the perturbation theory for the simplest case of the adiabatic approximation and 
apply them in section 6 for analysis of the evolution of the soliton parameters. A unique possibility is revealed to 
suppress the perturbation-induced shift of the soliton centre at the cost of proper matching of the soliton width and 
nonlinearity dispersion parameter. In section 7 wc calculate the spectral density of the radiation power emitted by 
the MNLS DM-soliton. Finally, section 8 contains a summary of the paper. 



II. BASIC EQUATIONS 

We model dynamics of a subpicosecond optical pulse in a cascaded transmission system with periodically varying 
chromatic dispersion by the MNLS equation: 

iq z + ^D(z)q„ + ia(\q\ 2 q) t + \q\ 2 q = iGq, (2.1) 



A' 



G = -7 + (cP* - 1) 5 ( z - n£). 

71=1 

Here all quantities are dimcnsionless: q{z,t) is the complex envelope of the electric field E normalized to the typical 
pulse power \E\ 2 = Po\q\ 2 , t is the retarded time normalized to the typical pulse width, the distance z down the fibre is 
normalized to the nonlinear length L n \ = (aPo) , a is the nonlinear coefficient. D(z) is the chromatic dispersion and 
a measures the strength of the nonlinearity dispersion. The loss coefficient 7 and a are supposed to be z-independent 
and I is the amplifier spacing. Taking q(z, t) in the form pEj 



we transform equation l|2.1|l to 



q(z, t) = u(z, t) exp 



dz'G(z') 



-D(z)u tt + iag{z)(\u\ 2 u) t + g(z)\u\ 2 u = 0, 



(2.2) 



where 



g(z) = exp 



2 ( dz'G(z' 
Jo 



As usually, we consider a transmission line comprising periodically alternating fibre sections and lumped amplifiers. 
Each fibre section of the length £ contains a portion of a compensating fibre with normal dispersion Z?_ < and 
length and a portion of a fibre with anomalous dispersion D + > 0. Amplifiers compensate for the fibre losses 
over the distance t. Hence, the total chromatic dispersion D(z) splits into two components: a small constant path- 
average (residual) dispersion S and rapidly varying local dispersion A(z). The local dispersion is periodic with 
zero average (A) = l~ x j Q dzA(z) = 0. The residual dispersion within the single fibre period is given by D lcs = 
[D-l- + D + (l — . We take the amplifier spacing £ to be equal to the DM map length. 

Different dispersion scales can be associated with the above dispersion map scheme: the dispersion lengths L± 
corresponding to the local chromatic dispersions D± and the dispersion length L rcs corresponding to the residual 
dispersion. We will consider the case of L + <C L Ies and L Ies ~ L n \ which corresponds to the condition £ -C 1. 
In accordance with this condition, it is natural to distinguish fast and slow scales which are represented by the 
variables £ = zjl and z, respectively. Therefore, 

D(z) = 5 + - e A(C). 

Following Ablowitz and Biondini we seek for a solution u(Cj z > t) of equation (|2.2|l as a series in powers of I: 



u(C, 2, t) = u (0) (C, 2, t) + fo (1) (C, z, t) + 



(2.3) 
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Writing the derivative d z as d z = £ 1 d{ + d z + 0(£), substituting the series l|2.3[l in equation l|2.2|) and equating terms 
with equal powers of i, we obtain a chain of equations. In the leading order we arrive at a linear equation 



,(0) 



(2.4) 



In the frequency domain (here and below we do not explicitly indicate the integration limits if they are infinite) 

u (0) (C,M) = J dLue- iuJt u (0 \(,z,Lu), u(°\C,z,u) = ^ / dte iwt w (0) (C,z,t) 
equation 12.4|l is easily solved: 



i / dC'A(C') + Co 



= U(°\z^)P(C,u). 



(2.5) 



Here Co = const. Hence, in the leading order the pulse evolution is mainly determined by the local dispersion A(^), 
while the integration constant £A°) (z, uj) depends on the slow variable z and is determined by the next-order equation: 



iuf + -A(C)u 



(i) 
tt 



In the frequency domain, where , z , t) = J duj exp(— iujt)u^ (£, z, uj), equation l|2.6|l has the form 



-A(C)fl« = -P(C,a;) 



W^iz^-^U^iz^) 



(2.6) 



(2.7) 



+g{z){l + au)) J J duj 1 duj 2 U {0 \z 7 cj + uj 1 )U {a \z 7 iv + Lj 2 )U {a) {z 7 uj + Lj 1 + uj 2 )T{(,uj 1 uj 2 ) 
where overbar means complex conjugation and 

T(C, k) = exp 



C 



i/t / dC'A(C') + C 



The standard procedure to remove secular terms, i.e., orthogonality of the right-hand side of equation l|2.7|l to a 
solution of the homogeneous equation, L d(P((,uj) • (r.h.s.) = 0, leads to the equation for U^(z, uj): 

{ Uf ] - -^-U {0) +g{z)(l + auj) / / duj 1 duj 2 U {0 \z,uj + uj 1 )U^\z,uj + uj 2 )U^ ) {z,uj + uj 1 + uj 2 )T(ujiuj 2 ) = 0, (2.8) 







T{k) = [ dCexp 


Hi 


JO 





(2.9) 



Equation (|2.8|l describes the averaged dynamics of the pulse with all fast and large variations removed and represents 
a straightforward generalization of the Gabitov-Turitsyn equation Q when accounting for the nonlinearity dispersion. 
It is easy to explicitly calculate the function T(k) (|2.9|l for the piecewise two-step dispersion map (Figure 0), given 

by 



A(C) = { A / ( 1 V 
CG "2-2 



s l 



2' 2 



where s = — A_(A + — A_) 1 is determined from the condition (A) =0. As a result, we obtain 

sin(^wicj 2 ) 



T{uj\UJ 2 ) 



\1UJ\UJ 2 
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FIG. 1: Two-step dispersion map with the residual dispersion 8 and local dispersions A±. 



Here \i = (1/4) [A+s — A_(l — s)] measures the normalized map strength. 

In what follows we will consider the case of small fi that corresponds to the limit of weak and perhaps moderate 
dispersion management. It gives 

T(wiwa) = 1 - -z^uful [i 2 <C 1. (2.10) 

D 

As a result, equation l|2.8|) takes the form 

iUc° } - -^-U {0) + g(z)(l + au) ( ( dwidw 2 L> (0) (z ,lu + uii + u>2) (2-11) 



= i/i 2 .g(z)(l + auj) J J &LUi&uj 2 U {a) {z,u + wi)L> (0) (z, u + w 2 )C/ (0) [z, uj + u\ + uj 2 )- 

Before we proceed further and come back to the temporal domain, we scale the coordinate z and the function U^: 
z = z'/S, U = {S/gf^U'. Performing the inverse Fourier transform and omitting primes, we arrive at the perturbed 
MNLS equation: 

iu z + iuf t + ia(\u\ 2 u) t + \u\ 2 u = ^ 2 (r + iar t ) (2.12) 

with the perturbation 

r(z, t) = u 2 u ttu + iuu t Uttt + 2(uuu + 2u 2 )u u + iu t u t tUt + Utt^- (2-13) 

It should be noted that the form (|2.13|l of the perturbation for the standard NLS-based DM soliton has been 
previously obtained by Lakoba and Pelinovsky [l4| and Boscolo et al. (24|. A contribution of the nonlinearity 
dispersion manifests itself by the a-dependent terms both in the left and right hand sides of equation (|2.12|) . In the 
next two Sections we will derive soliton solution of the MNLS equation and develop a perturbation theory to apply 
it to equation (|2.12|l . 

III. INTEGRABLE MNLS EQUATION 

The unperturbed MNLS equation 

™ z + ^Ww + ia{\u\ 2 u) t + \u\ 2 u = (3.1) 
admits the Lax representation A z — B t + [A, B] = with the matrices A and B of the form 

A = A(fc) 3 + 2ifcQ, Q = ( ° ) ' (3 - 2) 
B = n(k)(T 3 - 2fcA(fc)Q + 2ifc 2 Q 2 cr 3 + fcQ t er 3 - 2iakQ 3 . 
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Here 

A(fc) = ~h [Ak2 1); m = iA2(k) = _ i (4fc2 ~ 1)2j (3 - 3) 

k is a spectral parameter and 173 is the Pauli matrix. In the limit a — > equation i|3.1|) goes to the NL S eq uation. As 
regards the matrices A and B, the NLS limit is less trivial and is given by representing k in the form [26| 

k = i(l + aA) + C(a 2 ). (3.4) 

Then the limit a — > transforms the matrices A and B into the standard Lax matrices of the NLS equation with the 
spectral parameter A. 

A. Jost solutions 

The associated linear problem for the MNLS equation is written as 

J t = A(fc)[cr 3 ,J] +2ifcQJ (3.5) 

for the matrix-valued function J(t, k) (we will omit the evolution variable z, unless evolution aspects are concerned). 
Solutions of the spectral equation l|3.5[l obey two important symmetries, namely, the involution 

Jt(f,jfe) = J-^fc) (3.6) 

with \ denoting the Hermitian conjugation, and the parity 

3(t,k)=a 3 J(t,-k)a 3 . (3.7) 



Matrix Jost solutions J± of equation l|3.5[) obey the asymptotic conditions J± — ► 1 as t — > ±00. Since trA = 0, we 
have det J± = 1 for all t. Being solutions of the first-order equation 1)3.5(1 . the Jost functions are not independent. 
Indeed, they are interconnected by the scattering matrix S: 

J_=J+ESE- 1 , S(fc)=^~j)' detS = l, E = exp [A(fc)tcr 3 ] . (3.8) 

For the subsequent analysis, analytic properties of the Jost solutions are of primary importance. 

B. Analytic solutions 

Rewriting the spectral equation (|3.5|l for the Jost functions with the corresponding boundary conditions in the form 
of the Volterra integral equations, we obtain for the first column of the Jost matrix J_: 

k) (3.9) 



J- n (t,k) = l + 2ifc ( &b'u(t')J-2i{t', 

J — OO 

J_2i(£,fc) = 2ifc I dt'u(t',k)J- n (t',k)eKp 

J —OO 



-{Ak 2 - l)(t-t') 



It follows from the integrand in (|3.9I) that the column J_ is analytic in the domain C+ = {k, Imk 2 > 0}, i.e., in the 
first and third quadrants of the complex fc-plane, and continuous for Imk 2 = 0, i.e., on the real and imaginary axes. 
The same result follows for the second column 3 , of the Jost matrix J + . Therefore, the matrix function 



* + = (j W , J? 1 ) (3.10) 



solves the spectral problem (|3.5|l and analytic as a whole in C+. The analytic function $ + can be expressed in terms 
of the Jost functions and some elements of the scattering matrix: 

$+ = J+ES+E" 1 = J_ES_E-\ (3.11) 
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s + =(:?)- s += ss - 

Equation 13.11(1 yields 

det$+(t,fc) =a(k). (3.12) 
The function $+ can be expanded in an asymptotic series in A; -1 : 

$+(t,fe) =$f (tj + r^^fJ+O^ 2 ). (3.13) 
In virtue of the parity property l|3.7[) . expansion coefficients possess a definite structure. Namely, matrices <f>+ n \t), 



n 



— 0,1,,..., are diagonal, while 3>^ n+1 \t) are off-diagonal. Substituting the expansion (|3.13|) into the spectral 



equation 13.5fl yields the formula for the reconstruction of the potential Q: 



Q = l, 3 $«($f)- 1 , (3.14) 



as well as the equation for the leading-order term : 



$f t =-ia* 3 Q 2 $f. (3.15) 
Similarly we introduce a matrix function $Z (f, fc) made up from the rows (J^- 1 )^ of the matrices Jjj 1 : 

( J = X )m ^ 



^i 1 = 



2] 



This matrix is a solution of the adjoint spectral problem and analytic in the domain C_ = {fc, Imk 2 < 0}, i.e., in the 
second and fourth quadrants of the fc-plane, and continuous for Im fc 2 = 0. Besides, 





i _ 


ER+E ^J^ 1 




ER 






R_ = f 


1 




(! 


n 


b 


:) 



(3.16) 



R+= ( '.: , ] . R_ = I ; ' ! . R + S = R. 

Hence, det$I 1 (i, fc) = a(k) and the involution 

<f>Z\t,k) = & + (t,k) (3.17) 

takes place as well. 

C. The Riemann— Hilbert problem 

It follows from equation (|3.14|) that it is the analytic function $ + that determines a solution of the MNLS equation. 
On the other hand, the availability of two matrix functions and which are analytic in the complementary 
domains C± of the fc-plane and continuous on the common contour Im fc 2 = permits to pose the matrix RH problem 
for these functions: 

$I 1 (<,fc)$+(<,fc) = EG(fc)E~ 1 , (3.18) 



G(fc) = R+S+ = f J J V k e Im k 2 = 0. (3.19) 

Equation l|3.18|) is easily obtained from equations (|3.11|) and l|3.16|l . In other words, the RH problem l|3.18|l poses a 
problem of analytic factorisation of the non-degenerate matrix function G(fc) specified on the contour Im fc 2 = 0, in 
a product of two matrix functions (f> + and $I X analytic in the corresponding domains of the fc-plane. Hence, solving 
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the MNLS equation is equivalent to solution of the RH problem l|3.18[l . The normalization condition necessary to 
provide the uniqueness of solution of the RH problem follows from equation l|3.13[) : 

$+(t,k) -> $^ 0) (fc) at fc->oo. (3.20) 

Hence, we obtain the RH problem with the non-standard normalization ($+^ instead of the unit matrix), as distinct 
from the NLS case. 

In general, the matrices and $I X can have zeros kj and kg in the corresponding analyticity domains: 
det$>+(kj) = 0, kj G C+ and det$I 1 (fc^) = 0, kg £ C_. In virtue of the parity property l|3.7[l zeros appear in 
pairs, ±kj and ikg. The involution property l|3.17|) guarantees equal number J\f of zeros which we consider as simple. 
They are zeros of the RH problem that determine soliton solutions of the MNLS equation. 



D. Regularization of the RH problem 



We will solve the RH problem with zeros by means of its regularization, i.e., by extracting rational factors from 
<fr± that are responsible for the appearance of zeros. Namely, let us multiply $+(£, k) having a single zero kj, by a 
rational function k) which has a simple pole in kj. Then the product k)Sj (t, k) will be regular in kj. 

Since det = 0, there exists an eigenvector \xj) with zero eigenvalue, $+(kj)\xj) = 0. We take S" 1 in the form 

ET\t,k) = l + ^^Pj(t), P iW = MM, ( Xi | = | Xi > + , (3-21) 

where Pj is a projector of rank 1, = Pj. It is easy to see that detS" 1 = (kj — kj)(k — kj)~ x , as should be to have 
a pole in kj. The regularization in the point — kj is evidently given by 



where, due to the parity property (|3.7|) . \x~j) = a 3\Xj) an d P-j = ^^Pj^s- As a result, the matrix function 
E~ • is regular in ±kj. 



$+3,; E • is regular in ±fcj. Similarly, the matrix 2_^S^ $: x has no zeros in ±kg, where 



Regularising all 4M zeros of the RH problem, we represent the matrices $± as 

$± = 0±r, r = s_jv5v-S-iSi. (3.23) 

Here the matrices 4>±(t, k) solve the regular RH problem, i.e., without zeros, 

<jr l 4>+ = fege^f -1 , (3.24) 

while the rational matrix T(t, k) accumulates all zeros of the RH problem. Evidently, det T = ]Jj =1 (k - k 3 )(k - kj)- 1 . 
It follows from the structure of the rational multipliers S that the asymptotic expansion for T starts with the unit 
matrix: 

T(t,k) = i + k- 1 r (1) (t) + o(k- 2 ). (3.25) 

Besides, the matrix T hereditates the involution property T^(t, k) — r _1 (i, k). 

For the unperturbed MNLS equation, we can take solutions <p± of the regular RH problem (|3.24ll being k- 
independent. Hence, comparing the asymptotic expansions (|3.13(l and 1|3.25J) . we obtain that the solution of the 
regular RH problem is given by the leading-order term $ v} : 

<M*)=$f (t). (3.26) 
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In practice, a successive application of the elementary multipliers to regularize the RH problem would not be an 
optimal way. It is more convenient to decompose the product of Ej-'s l|3.23[) into simple fractions. After some algebra 
[23, H| we obtain 

T(t,k) = 1- £ (^SnlW^'^H (3-27) 

m , n — 1 

r-^t,*) = 1+ E ( fc -'OM( £, ~ 1 ) m ,>l 

m.n— 1 

with new vectors \m) (in the same way as |x±), the vectors |m) still solve the equation $ + (K m )|m) = 0) and matrix 
D with the entries 

(n I m) 

D nm = V ^1 . (3.28) 

We adopt the following enumeration of zeros: K m = {k%, —hi, k 2l —k 2 , ■ . . , fcyv, — fcjv}. 

To determine coordinate dependence of \m), we turn to the equation $ + («; m )|m) = 0. Differentiating it in t and 
accounting equation (|3.5() yields |m)t = A(K m )<73|m). In the same way, differentiation in z together with the evolution 
equation 

= B$+ - 0(fc)a 3 $+ (3.29) 

gives \m) z = Q(K>m)c3\m) ■ Therefore, the coordinate dependence of the vector |m) has the very simple form: 

\m) = exp{[A(K m )t + Q(K m )z] cr 3 } |m ), (3.30) 

where \mo) = (pi,P2) T is a constant vector. The parity property (j3.7(l gives |2m) = cr3|2m — 1). 

Zeros K m and vectors \m) comprise the discrete data of the RH problem that determine the soliton content of a 
solution to the MNLS equation. The continuous RH data is characterized by the matrix G(k), k € Im k 2 (or by the 
scalar function b(k) entering G(fc) (JHZE50 an d is responsible for the description of radiation. The evolution equation 
for G(fc), 

G z =fi(fc)[o- 3 , G], (3.31) 

easily follows from equations l|3.18|l and l|3.29|l . 

E. Reconstruction of the potential 

The rational function T(t, k) plays the key role in both finding solitons and describing its perturbed dynamics in 
the adiabatic approximation. Indeed, equations H3.5(l . Ij3.14|l and (|3.23(l give 



iaa 3 Q 2 Uf) 1 $+ + « 0) ) 1 (A(fc)K *+] + 2ifcQ$+) 



r t - 

Hence, the equation for T^ 1 = [T(k = 0)] _1 , 

(r - 1 ) ( =-ia<T 3 Q 2 To\ (3.32) 

coincides with equation l|3.15|l for <I>^ ' ) . Identifying both functions, 

ro 1 = $<?>, (3.33) 
we derive from Ij3.14|l and (|3.26(l a formula for soliton solutions in terms of T: 

q = |o- 3 ro x rWr . (3.34) 

Here T^ 1 ^ is the coefficient in the series (|3.25f) . Note that equation (|3.34|) cannot be applied to consider radiation 
emitted by a perturbed soliton. 
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F. MNLS soliton 



Now we apply the formalism developed in the preceding Section to get a single soliton solution of MNLS. The RH 
data are purely discrete: TV = 1, b(k) = (or G(fc) = 1), + = <f>_. As it follows from (|3.34|) . to obtain the soliton, 
we should know the rational function Y which in turn is completely determined by zeros K\ and k 2 (i.e., by ±fci), as 
well as by the vectors |1) and |2) = 03 11). From ll.'l-'ill and (|3.30|l we obtain the vector |1) explicitly: 



1\ — P (l/2 (to+iVo f e ] 
^ _ e I -1/2 r+i^ I 



Here r and ■0 are natural soliton coordinates related linearly to i and z: 

T 

1> 



w(t -t), 

V ~ 

-0, 



w 



t = VZ — To/w, 

%/} = -(v 2 + w 2 )z + ip , 



(3.35) 



(3.36) 



"00 



-T + 1p - 



The real parameters v and w are constants which are determined by the zeros k% and ki and will be identified with 
the soliton velocity and inverse width, respectively: 



v = -\i- 2(kj + m 

a 



u; = -(*?-*?); 



(3.37) 



To and ipo are real parameters which are given by the components of the constant vector |1q) (|3.3Ufl . p\/p2 — exp(ro + 
i^o)- Note for the future use that the constant soliton parameters acquire in general a slow ^-dependence in the 
presence of perturbation. This results in a modification of equations 



t : 



1 

w 



dz'v(z')w(z') — tq ] , ip = vi 



dz' (v 2 (z')-w 2 {z'))+ij a . 



Having explicitly the vector |1), we find from l|3.27|l the matrix T 

D D+ 



r = 1 



k — k\ k + k\ 



v- 1 = 1 



D 



k — k\ k + k\ 



(3.38) 



(3.39) 



where 



D 



k\ - k\ 



c+ 



D 



k\ - k\ 



C+ 



C+ c+ 



D 



-cr 3 D_(T 3 , 



-(T 3 D-<T 3 , 



and the functions c± are given by 

c± = fcic ±r + k x e TT . 
Substituting this Y in (|3.34|l gives a simple expression for the MNLS soliton: 



u s = — -5-e 

1 ci 



i-0 



(3.40) 



(3.41) 



This soliton solution, like that of the NLS equation, depends on four real parameters: velocity v, inverse width w, 
initial position To and initial phase i/jq. Nevertheless, its properties essentially differ from those of the NLS soliton. 
Indeed, the soliton l|3.41() has no habitual sech-like shape though its envelope is close to the hyperbolic secant. The 
square of module has the form 

1 



M 2 = ^ 2 



1 — av + — av) 2 + a 2 w 2 cosh(2u>(T — vz)) 
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Hence, the envelope \u s \ moves holding its shape and there is no any self-steepening, contrary to the behavior of the 
NLS soliton under the action of the nonlinearity dispersion. The reason lies in the phase properties of the MNLS 
soliton. Namely, we can write u s as \u s \ exp(h/> + iip n i), where 

( n i 3 — tan 2 6 tanh 2 r \ , ifl 

ip n l = arctan tan tanh r * — , fei = |fci|e . 

V 1 — 3 tan 2 6 tanh r / 

We see that the MNLS soliton is characterized by a highly nonlinear phase and hence is intrinsically chirped. It is 
this nonlinear phase that prevents the MNLS soliton from distortion of its shape, as opposite to the NLS soliton 
whose linear phase cannot withstand the self-steepening. Further, because w ~ a -1 , the soliton (|3.41(l contains the 
nonlinearity dispersion parameter a in the denominator. Thereby, the soliton (|3.41|) is non-perturbative in a and 
cannot be obtained in the framework of the NLS equation with a-dependent perturbation. At the same time, in 
accordance with the limit procedure l|3.4[) . the MNLS soliton reproduces the standard NLS soliton in the limit a — > 0. 
Besides, the optical energy 

£= [ dr\u s \ 2 = -6 (3.42) 
J a 

being the invariant of the MNLS equation, has the upper limit 2ir/a because 6 < ir/2. 

IV. PERTURBATION THEORY FOR THE MNLS SOLITON 

In this Section we perform a general analysis of the perturbed MNLS equation 

iu z + -u tt + ia{\u\ 2 u) t + \u\ 2 u = eR. (4.1) 

Here R determines a functional form of a perturbation (e.g., equation l |2.13| l). e is a small parameter. In general, a 
perturbation forces the RH data to slowly evolve in z. Indeed, a perturbation causes a variation <5Q of the potential 
entering the spectral equation H3.5|l and hence a variation of the Jost solutions: 

6J± t = A(fc)[a 3 , SJ±] + 2ifc((5QJ± + Q6J±). (4.2) 

Solving H4.2H gives 



S3 ± = 2ifcJ±E ^ dt'F,- 1 J ± 1 dQJ±~Ej E _1 . 



(4.3) 



To distinguish between the 'integrable' and 'perturbative' contributions, we will assign the variational derivative 5/5t 
to the latter. Then, evidently 

^ = eR ' R= ( -R R o 

By means of equations l|3.8|l , (|3.11|) and (|4.3|) we find a variation of the scattering matrix: 

^ = 2e/cS + (^J dtE-^^Rfc+EJ Si 1 = 2ekH+ 1 (J dtE- 1 $I 1 R<I>_E^ R . 
It should be stressed that they are the analytic solutions <I>± that naturally enter this equation. Let us denote 

nt-2 

Y±(h,t 2 ) = 2k i diE^S^RSiE, T±(fc) = T±(-oo, oo). (4.4) 
Jti 



Then 



— = eS+T + (fc)Sl 1 = eR7 1 T_(fc)R_, (4.5) 

oz 

where the matrices S± and Rj- have been defined in equations (|3.11() and l|3.16|l . The matrices T± are interrelated 
by the matrix G entering the RH problem Ij3.18|l : 

T_(fc) = GT+G" 1 . (4.6) 
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Finally, the variations of the analytic functions <f>± follow from (|3.25|) and i|3.30[l : 



5<P, 



e$ + En + E _1 , 



5$: 



Sz Sz 
Here the evolution functionals H±(t, k) are defined in terms of T± p9|: 



eEn_E _i $ 



U+(t,k) 



T+n(— oo, t) -T + i 2 (t, oo) 
T +2 i(-oo,t) -T +22 (i, oo) 



n_(t,fc) 



T_u(-oo,f) T_i 2 (-oo,t) 
-T_2i(t,oo) -T_ 22 (i, oo) 



(4.7) 



The evolution functionals contain all needed information about a perturbation and enter the evolution equations for 
the perturbed RH data (see below). As a result, the evolution equations for $± gain additional terms responsible for 
the perturbation: 



4> 



+ z 



B$+ - <f> + n(k)<r 3 + e^+En+E" 1 , 



(4.8) 



®zl = -$: 1 B + r2(fc)cr 3 $z 1 -eEn_E" 1 $: 1 . 

Then from equations l|3.18[) and (|4.8I) we obtain the evolution equation for the perturbed continuous RH data: 



G z = £l(fc)[er 3 , G] + e(Gn + — n_G), Imfc 2 



0. 



(4.9) 



Since the left-hand side of 1)4. 9|) does not depend on t, we can consider this equation for t — > oo. This gives simplified 
formulas for the evolution functionals: 



n+(fc) 



T+n(fc) 0' 

T +2 i(fc) or 



n-(fe) 



T-n(fc) T_ 12 (fc) 




(4.10) 



To derive evolution equations for the discrete RH data, let us turn to the equation $ + (fci)|l) = which is valid 
irrespectively of the presence of a perturbation. Taking the total derivative in z, we get 



az 



|i> + Mfci)^|i) = o. 



Assuming a perturbation-induced z-dependence of k% yields 

^$+(fc) = B$+ - fl(k)^+a 3 + e^+En+E- 1 + k z ^+(k). 

Recall that the evolution functional n + (fc) is expressed in terms of elements of the matrix T + (see l[4.7jl ) which in 
turn depends on 4>^ 1 (see l|4.4|l '). Therefore, T+ has the simple pole in k\ and hence n + (fc) represents a meromorphic 
function with simple poles in zeros of $ + . This implies that near k\ we can write 



n + (fc) = < og >(fc) 



k — k\ 



Res[n+(fc), kx] 



where n!f g) (fc) is a holomorphic part of n + in the point k\ . Following then the reasoning of 30] , we obtain that the 
perturbed evolution of the vector |1) is given by 

\l) z - n(fci)<T 3 |l) - eEikx^^ikx)^!}. 
In terms of the i-independent vector |1) = E _1 (fci)|l) the above equation takes the form 



|i), = n(fci)o- 3 |i) 



?txf(kx) 

rtf(kx) o 



|i>- 



(4.11) 



As regards the evolution equation for kx, it follows by taking the total z-derivative of det = and has the 

form 



kxz 



d z (det$+(fci)) 



9 fe (det 



|fci 
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It has been shown (see l|3.23|l ) that $ + = <j>+T. Because V — S_iSi, we have detT = (k 2 — kf)(k 2 — kf) Besides, 
equations l|3.26|) and 13.33fl gives cf) + = T^ 1 (k = 0) and det(j> + = k\/k\. Hence, 

Taking into account the relation d z [det $ + (/«)] = e (trll + (k)) det $+(£;) and equation l|4.10|l . we eventually arrive at 
a simple evolution equation for the zero ki : 

k lz = -eRcs[trn + (fc),fci] = -eRes [T +n (k), h] . (4.12) 

Hence, equations H4.9J) . (|4. 1 1|> and 1(4. 12j) determine perturbation-induced evolution of the RH data. These equations, 
however, cannot be directly applied because n-t- and T± entering them depend on unknown solutions <f>± of the spectral 
problem with the perturbed potential Q. The smallness of e allows us to develop the iterative scheme to consecutively 
account for two basic approximations: the leading-order adiabatic approximation and the next-order one. 



V. ADIABATIC APPROXIMATION 



In the framework of the adiabatic approximation we assume that the perturbed soliton does not radiate and adjusts 
its shape to the unperturbed one at the cost of slow evolution of its parameters. Hence, only the discrete RH data 
are relevant in this approximation. The soliton coordinates are given by equations l|3.36[) and (|3.38() . It is seen from 
(14.11(1 and l(4.12|l that the matrix T+ completely determines evolution of the parameters. Because now = r o " 1 r, 
we obtain in the adiabatic approximation 



T+(fc) =2k J diE^r^roRT^rE. 
The explicit form of T is given in l|3.39|l . Therefore, 

Hes[T +u (fc),fe 1 ] = -iafc? / £ (p(r) + p(-r)) e 



(5.1) 



where p(r) = Re ^ . Then equation l|4.12|l gives a very simple equation 

kx z = -eak 1 / -j-p+e , 



(5.2) 



where p± = p(r) ± p(—r). In addition, we can derive from (|3.37|) evolution equations for the soliton velocity and 
inverse width: 



-2ie 



dr 



(fc x V-fc?e- T ) p+(r), 



w, = 2e 



At 



(5.3) 



Somewhat more cumbersome calculation is needed to derive the evolution equation for the vector |1). In accordance 



with [j4~TT)l we need know TY eg) (/ui). Hence (R = r RIV 1 ), 



T+(fc) - 



Res [T + (fc),fci] 



k — k\ 



1 



= 2Ai / dtE- 1 ^) I 1 + — D+ ) HT(fci)E(*i) 

\k— k± J \ 1 



8i r f 
+—kf cr 3 , / dttE" 1 (fci)D_Rr(fc 1 )E(fci) 
a [ J 

Taking now into account that 

1) = exp 



2/ dfE _1 (fci)D_R 1 



fc x D_ 



(fci-fci) 2 (h + ^y 



E(fci). 



dz'fi(fci)<7 3 



l(Pi,Pa) r > 



13 



and P1/P2 = exp(ro + iipo)-> we nn d from (|4.11() 



TOz 



dr 



-ac + - — (kle T + kfe~ T ) 



dr 



w 



^ z = 2iet / — (k\e T - kfe- T ) p+ - i 



I I dr 



2t (k\e T + kle~ T ) 



(k 2 + kf) c+ + |/ci| 2 c_ - 2r (kle T - k\e 



P-- 



Finally, differentiation of t and ip (|3.38[) in z gives equations for the soliton centre t and phase 

e f dr 



v-2- 



-ac+ (k?e T + fc?e~ T ) 

s w v ' 



P- 



tpz = ~ (v 2 + w 2 ) + — I %■ I — [kl(v + iw)e T + k\{v - iw)e T ] - - i\k r \ 



C_ > p 



(5.4) 
(5.5) 



It is remarkable that the even perturbation (p_ = 0) influences the soliton velocity and width only, while the odd 
perturbation (p + = 0) does not modify these parameters and does modify soliton position and phase. 

VI. SUBPICOSECOND DM-SOLITON IN THE ADIABATIC APPROXIMATION 



Now we invoke equations l|5.3[l - l|5.5[l to analyze the DM-soliton parameters in the adiabatic approximation with 
e = p 2 and R = r + iarj. Substituting the soliton solution u s l|3.41[l into 12. 1311 . we obtain 



r = -16iw 7 |fc 1 | 4 4^e i ^, 



(6.1) 



where 



f(r) 



6 ( 4 
- + 12 



Jlr i „-2t 



1 f 1 



(6.2) 



c+c_ \ ci c 

It follows immediately from (|6.2() that 

p+ = [{r + ior t )(r) + (f — iaf t )(— r)] exp(— iip) = 0. 

Hence, the DM-soliton preserves its velocity v and width w^ 1 in the adiabatic approximation. Without loss of 
generality, we pose in what follows v = 0. 
As regards p_, it is represented in the form 



P- 



-64iw 7 |fci| 4 fc 1 



(l+/9 2 )i-(l-/3 2 )fi 



where (3 — k\jk\ = exp(2i0). In fact, we need not calculate the term (//c3_) T because the integration in r in l|5.4|) 
permits to transform the integral with f Tl by means of integration by parts, into integrals with /. As a result, we 
obtain the following formula for t z : 



t z = -32^aw b |fc 1 | 4 /fc 4 



dr 



{l + 4/? 2 + /3 4 + 3/3 2 (f3e 2T + /T^" 2 ") 



(6.3) 



+ - 



4r 



1-/3 2 

Calculation of integrals in l|6.3|) is straightforward and gives 

t z = S2p, a 



-1 + 2f3 2 + 2/3 4 -f3 6 + f3 3 (f3 2 e 2T + /r 2 e~ 2T )] }/. 



o 



' P l( 0)-i^P 2 (0) 



12 Jfei 



aw 



(6.4) 



Pi (9) = 7 cos(8<9) + 244 cos(46>) +278, P 2 (0) = 8 cos(8<9) + 35 cos(46>) + 44. 



(6.5) 



14 



0.5 



theta / 
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FIG. 2: The curve 16.71 intersects the horizontal axis in the point 8\ — 0.138. The soliton width and the nonlinearity dispersion 
parameter a obeying the condition aw — tan(2#i) provide the perturbation- free- like propagation of the DM-soliton. 



The condition v = implies in accordance with 1|3.37(1 and k\ — \ki \ exp(i(9) the following expressions for |fci| 2 and w: 

4|A:i| 2 = {l + a 2 w 2 ) 1/2 , aw = tan(20). 
Therefore, equation (|6.4|l takes the form 

Now it is reasonable to take up a question of how to minimize a perturbation-induced deviation of the soliton centre 
position from that of the unperturbed soliton. Evidently, t z will be strictly zero in zeros of the curve 

Y(6) = 249P 2 (6) - Pi {6) sin(40). (6.7) 

The plot of the curve 1(6.7(1 is shown in Figure El It is seen that there exists indeed the nontrivial point 9\ of 
intersection of the curve with the horizontal axis. Numerical solution of equation ((6.7(1 gives 0\ — 0.138. It means 
that if the nonlinearity dispersion a and the inverse soliton width w obey the condition aw = tan(2#i) = 0.273, the 
perturbation would not disturb dynamics of the soliton centre. Note that such a compensation is impossible for the 
NLS DM-soliton. 

Just in the same manner we can calculate a soliton phase variation from equation 1)5.5(1 . However, for the single 
soliton this parameter is less informative as compared with the soliton centre variation, and we do not exemplify here 
the results. 



VII. PERTURBATION-INDUCED RADIATION BY DM-SOLITON 



For treating radiation effect, we should go beyond the adiabatic approximation. It means that we should abandon 
the condition G = 1 valid in the adiabatic approximation and find the function b(k) entering the matrix G l(3.19|l . 
Indeed, the spectral density of the radiation energy emitted by the perturbed soliton is given by p3| 

Wfc) = -^7^1n[l-sgn(fc 2 )|fe(fc)| 2 ] • 
The function b(k) accounting for radiation is of the order of /j, 2 . Hence, for /i 2 C 1 we have 

Then the spectral density of radiation power is written as 

W(k) = ^-£ iad (k) = — ^Re(6 b z ). (7.1) 
az ira \k\ 
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FIG. 3: Typical dependence of the spectral density of radiation power W(k) (in arbitrary units) on wave number. Two peaks 
correspond to forward and backward radiation emitted by perturbed soliton. 



Now we turn to the evolution equation (|4.9|l for the matrix G which determines the RH problem (|3.18|1 . Taking 
the evolution functionals IT± in the form (|4. 10ft . we immediately obtain the equation for b: 



b z = 2Q(k)b- M 2 T+i 2 . 



(7.2) 



Note that T + = Y_ in the order of /i 2 , in accordance with equation l|4.6(l . The matrix element T + i2 is calculated via 
(|5.1|) where we drop the terms which rapidly decrease at infinities: 



, k 3 (l + a 2 w 2 ) 3 e 

T+i2 = 4tt— = ^-e e cos 2(9 exp 

\k\ r a'w 



-(fc 2 -l)^+i0o 
2 w 



exp ( —w zj H(k,9)sech-^ 



where 



H(k, 9) = gki {2p [-17 - p 2 + (p 2 - 3) cos(40)] - 2i [-19 - 3p 2 + (3p 2 - 1) cos(46>)] } 
- g' x k x {p [59 - p 2 + {p 2 - 79) cos(40)] + 5i [9 - 3p 2 + (3p 2 - 29) cos(46»)] } , 



(7.3) 
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k 2 -k\ 
k 2 -k 2 ' 



4k 2 



Let us remind that aw = tan(2#). Inserting l|7.2|l and (|7.3I) into (|7.1I) . averaging in z over the chromatic dispersion 
map period I and dropping fast oscillating term, we obtain the following expression (of the order of C(^i 4 )) for the 
spectral density of radiation power: 



W(k) = - 
. a 



id 



TTfi 4 k 4 



(l + a 2 w 2 ) 4 



iw 2 [a 2 w 2 + (4k 2 - l) 2 ]' 



\H(k, 



|2„20 



e 2W sech" [ - 

2 aw 



7T 4fc 2 



(7.4) 



Figure demonstrate a typical dependence of the radiation power on the parameter k (which is proportional to the 
wave number). It has two peaks which correspond to forward and backward radiation. At the same time, radiation 
power depends only slightly on the nonlinearity dispersion parameter a, this dependence being detectable for a > 1. 



VIII. CONCLUSION 



In this paper, dynamics of a subpicosecond optical pulse in a dispersion-managed fibre links has been studied. 
The analysis has been performed in the limit, where the nonlinearity and average chromatic dispersion affect the 
pulse evolution on a much larger distance than the local dispersion does. We have derived the Gabitov-Turitsyn 
equation for the subpicosecond DM-soliton amplitude. For the case of the weak chromatic dispersion management, 
this equation has been reduced to a perturbed modified NLS equation. In the framework of the RH problem associated 
with the MNLS equation, a perturbation theory for the MNLS soliton has been elaborated. A possibility has been 
discovered to control the perturbation-induced shift of the soliton centre position choosing properly the soliton width 
and nonlinearity dispersion parameter. Spectral density of radiation power emitted by the DM-soliton has been found. 
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We believe the integrable MNLS equation serves as a true basic model to account for the Raman self-frequency 
shift and third-order chromatic dispersion. Both these effects have been analysed for the NLS DM-soliton in papers 
[sTL l8^ | . An example of the Raman effect action on the usual MNLS soliton was considered in . As regards the 
third-order dispersion effect for the MNLS soliton, the standard perturbation theory cannot capture essential features 
of the soliton dynamics, and a perturbation method that goes beyond all orders is needed (see [34j for the NLS 
soliton) . This challenging problem is outside the scope of the present paper. 

Acknowledgments 

This work was partly supported by the Belarussian Foundation for Fundamental Research under Grant No. $04P-071. 



[1] N.J. Smith, F.M. Knox, N.J. Doran, K.J. Blow and I. Bennion, Electron. Lett. 32, 54 (1996). 

[2] I.R. Gabitov and S.K. Turitsyn, Opt. Lett. 21 327 (1996); JETP Lett. 63, 327 (1996). 

[3] I. Gabitov, E.G. Shapiro and S.K. Turitsyn, Phys. Rev. E 55, 3624 (1997). 

[4] S.Kumar and A. Hasegawa, Opt. Lett. 22, 372 (1997). 

[5] T.I. Lakoba and D.J. Kaup, Electron. Lett. 34, 1124 (1998). 

[6] S.K. Turitsyn, T. Schafer, N.J. Doran, K.H. Spatschek and V.K. Mezentsev, Opt. Commun. 163, 122 (1999). 

[7] S.B. Medvedev, O.V. Shtyrina, S.L. Musher and M.P Fedoruk, Phys. Rev. E 66, 066607 (2002). 

[8] A.Labruyere, P. Tchofo Dinda and K. Nakkeeran, J. Mod. Opt. 52, 2731 (2005). 

[9] B.A. Malomed, Soliton Management in Periodic Systems (Springer, New York, 2006). 
[10] C. Lin, H. Kogelnik and L.G Cohen, Opt. Lett. 5, 476 (1980). 
[11] M.J. Ablowitz and G. Biondini, Opt. Lett. 23, 1668 (1998). 

[12] V.E. Zakharov, in Optical Solitons: Theoretical Challenge and Industrial Perspectives, edited by V.E. Zakharov and S. 

Wabnitz (Springer- Verlag, Berlin, 1999), pp. 73-90. 
[13] Yu.V. Lvov and I.R. Gabitov, Dispersion management in optical fiber links: Integrability in leading nonlinear order, 

arXiv:chao-dyn/9907007. 
[14] T.I. Lakoba and D.E. Pelinovsky, Chaos 10, 539 (2000). 

[15] Yu.S. Kivshar and GP. Agrawal, Optical Solitons: From Fibers to Photonic Crystals, (Academic, San Diego, 2003). 

[16] A. Rangwala and J.A. Rao, J. Math. Phys. 31, 1126 (1990). 

[17] E.V. Doktorov and V.S. Shchesnovich, J. Math. Phys. 36, 7009 (1995). 

[18] A. A. Zabolotskii, Zh. Eksp. Teor. Fiz. 120, 252 (2001). 

[19] Yu.S. Kivshar, Physica D 40, 11 (1989). 

[20] E.V. Doktorov and R.A. Vlasov, J. Mod. Opt. 38, 31 (1991). 
[21] V.S. Shchesnovich and E.V. Doktorov, Physica D 129, 115 (1999). 
[22] X. Chen and J. Yang, Phys. Rev. E 65, 066608 (2002). 
[23] V.M. Lashkin, Phys. Rev. E 69, 016611 (2004). 

[24] S. Boscolo, J.H.B. Nijhof and S.K. Turitsyn, Physica Scripta 62, 479 (2000). 

[25] A. Hasegawa and Y. Kodama, Opt. Lett. 15 1443 (1990); Phys. Rev. Lett. 66, 1611 (1991). 

[26] E.V. Doktorov and I.S. Kuten, Europhys. Lett. 53, 22 (2001). 

[27] T. Kawata, in Advances in Nonlinear Waves, edited by L. Debnath (Pitman, London, 1984), p. 210. 
[28] E.V. Doktorov and V.M. Rothos, Phys. Lett. A 314, 59 (2003). 
[29] V.S. Shchesnovich, J. Math. Phys. 43, 1460 (2002). 

[30] E.V. Doktorov, N.P. Matsuka and V.M. Rothos, Phys. Rev. E 68, 066610 (2003). 
[31] T.I. Lakoba and GP Agrawal, J. Opt. Soc. Am. B 16, 1332 (1999). 
[32] T.I. Lakoba and D.J. Kaup, Opt. Lett. 24, 808 (1999). 

[33] A. A. Afanas'ev, E.V. Doktorov, R.A. Vlasov and V.M. Volkov, Opt. Commun. 153, 83 (1998). 
[34] P.K.A. Wai, H.H. Chen and Y.C. Lee, Phys. Rev. A 41, 426 (1990). 



